Quantum phase coherence in non-Markovian and reaction-diffusive transport 
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We study quantum phase coherence and weak locahzation (WL) in disordered metals with re- 
stricted back-scattering and phenomenologically formulate a large class of unconventional transport 
mechanisms as modified diffusion processes not captured by the Boltzmann picture. Inspired by 
conductivity measurements in ferromagnetic films and semiconductors where anomalous power law 
corrections have been observed, we constrain memory dependent, self avoidance effects onto the 
quantum enhanced back-scattered trajectories, drastically altering the effect of weak localization in 
two dimensions (2D). Scale dependent corrections to the conductivity fail to localize the electrons in 
d > 2 for sufficiently weak disorder. Additionally, we analyze quantum transport in reaction-diffusion 
systems governed by the Fisher's equation and observe asymptotically similar delocalization in 2D. 
Such unconventional transport might be relevant to certain non-Fermi liquid or strongly correlated 
phases in 2D within the negative compressibility regime. 

PACS numbers: 73.20.Fz, 72.15.Rn, 71.10.Hf, 47.70.Nd 



Introduction- Lacking an adequate mathematical the- 
ory of the transition point and no rigorous proof of 
the existence of extended states in three dimensions 
(3D), Anderson localization of non-interacting electrons 
is an ongoing endeavor of condensed matter physics fH]. 
Modern approaches exploit symmetry classes of a given 
disordered Hamiltonian, constructing perturbative con- 
ductance g scaling forms of the diffusion (Goldstone) 
modes. 10 such universality classes have been identi- 
fied with certain cases demonstrating an absence of lo- 
calization in lower dimensions such as the symplectic 
class, which describes systems with extrinsic spin-orbit 
interactions 0, 13] • Weak localization (WL), is a powerful 
comphmentary approach, which in contrast to Ander- 
son or strong localization, are its precursor effects aris- 
ing from phase coherent back-scattering Q in the metallic 
regime. 

Following its early development, there has been much 
progress in understanding how the quantum phase co- 
herence is affected by fundamental processes such as 
extrinsic and intrinsic spin-orbit interactions, spin-flip 
scattering, electron- phonon interactions, magnetic fields, 
Nyquist noise etc....[4| These earlier cases demonstrated 
a renormalization of the diffusion constant, which can 
be realized in important phenomena such as anoma- 
lous magneto-conductance, anti-localization and spin- 
selective localization;2, iSi, (6|. However the existing theo- 
retical studies have not successfully addressed the conse- 
quences of non-equilibrium conditions, strong interaction 
ground states or non-Markovian processes. Therefore it 
is our purpose here to initiate an exploration of trans- 
port paradigms beyond the ordinary case of simple dif- 
fusion, by including certain observed features of strongly 
correlated groundstates such as possible non-equilibrium 
conditions and non-Markovian affects. 

Interactions have been included in disordered elec- 
tron systems both from a direct perturbative framework 
(Hartree Fock -I- impurity vertex corrections) f?! , a non- 



linear sigma model app roach[i-[i3] and from an alterna- 
tive scaling analysis [111]. More recently, the experimental 
discovery of a 2D metal-insulator transition (2DMIT) in 
high mobility SI-MOSFET hetero-structures has received 
much attention lately fl2l|. Following its discovery, a 2 
parameter scaling (3 function which takes the limit of a 
large valley degeneracy n.y — )■ cx) in a 1 /n^ expansion has 
yielded important and favorable results in understanding 
transport in dirty Fermi liquid metals[l3,[l4[. 

Experimental evidence suggests that the transition 
point of the 2DMIT takes place within a negative com- 
pressibility regime near Ts « 10 for a wide variety of elec- 
tron and hole carrier systems 15]. It should be noted that 
the neutralizing background renders the total compress- 
ibility positive, however this contribution is subtracted 
off and the density-density correlation function in the 
limit (7 ^ for the electron gas itself is negative. In 
addition, it is well known that strong correlation effects 
may lead to negative compressibility ground states, al- 
though the proper fixed point of the correlated fluid in 
2D has not been established. A. Schakel has shown that 
such conditions may lead to analogues of charge density 
wave formation and droplet clustering consistent with 
local density inhomogeneities or imbalance 16| . There- 
fore with these experimental and phenomenological con- 
ditions in mind, and as an alternative approach we for- 
mulate the problem as a reaction-diffusion diffusion pro- 
cess, predicting the destruction of phase coherence and a 
possible description of non-Fermi liquid groundstates in 
the negative compressibility regime, not captured by the 
ordinary Boltzmann picture. 

Another interesting case of unconventional transport 
involves non-Markovian effects, which is a manifesta- 
tion of correlated disorder. Earlier studies have detected 
memory dependent effects by applying an external mag- 
netic field to a 2D electron gas, with observable conse- 

As another application 



quences in the transport |17Hl 
of the semiclassical formalism we include infinite repul- 
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FIG. 1. Typical back-scattered trajectories and the expected 
time dependent behavior of their end to end distances R{t). 

(a) In a system containing self avoidance none of the paths 
taken by the electrons may be taken twice, causing the effec- 
tive scattering length to be dependent on the dimensionality. 

(b) Ordinary diffusion, where each scattering event is statis- 
tically independent of previous ones, allowing for the crossing 
or retracing of steps. 



sion of the time reversed interference trajectories which 
can be incorporated into the WL corrections through the 
scale dependence of the effective scattering length R{t). 
In simple terms the scattered electrons will remember 
their past and avoid taking trajectories already traversed, 
known as the extreme limit of repulsion within the ran- 
dom walk called self avoidance. 

Semi- das steal weak localization- WL can be under- 
stood semi-classically in terms of the diffusive behavior 
of a particle in a d-dimensional disordered system. Typi- 
cally, when the disorder is weak, the mean free path Z, is 
much greater than the quantum-mechanical wavelength 
Xp oc kp^, kpl » 1 such that the conductivity is modi- 
fied by. 



o- = (To + Sa, \5a\ <C ctq 
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with the usual Drude conductivity given by (Tq = 
e'^pr/me for charge e, carrier density p, scattering time 
T and mass mg. It is generally appreciated that WL is 
a quantum-interference contribution to 6cr that results 
from phase-coherent back-scattering. Qualitatively, the 
quantum interference corrections can be understood from 
within a Feynman path integral description 0, [IS] ■ The 
Feynman paths (indices i,j) connecting two locations 
within the disordered metal can be described by prob- 
ability amplitudes Ai. The total probability to reach the 
two points W must sum all possible weighted trajecto- 
ries, including a cross interference term W = ~ 
J2i + ^ back-scattered trajectory that 

interferes with its time reversed path is quantum mechan- 
ically enhanced and facilitates localization^]. 

As a practical matter, in order to calculate 5a we con- 
sider some variant of a diffusive transport along a closed 
trajectory, fixing the particle's classical return probabil- 



ityCD(f, 01201, 
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dt 
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where tpi^D, "r, t) is a function that depends on the spe- 
cific Hamiltonian considered. Note that C^i is expected 
to be a normalized probability distribution such that non- 
singular averages of length {x) , (x^) can be computed. 
Therefore there should be present a proper competition 
of terms (growth and dissipation) in ■0(CD,r, i). Fur- 
thermore, in systems with density inhomogeneities and 
subsequently nontrivial -(/'(Cd, r, i), it is not essential for 
maintaining a consistently normalized Cooperon that the 
fiow of matter in and out of a local space be equal. 

Upon determining CD{r,t) the 6a correction is ex- 
tracted from classical considerations of a compact vol- 
ume i?'*(over the back-scattered region), spreading as a 
function of time. R{t) is a central quantity in the study 
of random walks, which can accurately describe the dis- 
crete scattering events of the metallic regime. In the 
continuum limit, the random walk becomes a diffusive 
process controlled by CD{r,t). Therefore, R{t) is deter- 
mined from C'd {r, t) by the relation, 
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which is essentially the root mean squared value of the 
position taken by the electron. To determine the full con- 
tribution to (5(7, R{t) must be integrated in the following 
manner, 



2e2 
7:n 



dtR{ty 



(0.4) 



where the probability of return can be expected to be 
damped R{t)~'^ — R{t)~'^e~^*'/'^'^\ We therefore include 
\/t,^ as the dephasing rate that regularizes the integra- 
tion. The interesting case is given by ^ r. Note 
that the integral in Eq. (j0.4|) is precisely the result of the 
maximally crossed diagrams or Cooperon-like terms gen- 
erated in the formal Green's function treatment. For 
the ordinary case of non-interacting electrons in zero 
external fields, r^p ~ Tie for a given inelastic scatter- 
ing time Tie and ip{CD,r,t) = 0, yielding Coifjt) = 
g-|r-ro|V4Dt/(4^£)^)d/2 ig^diug to R{t) = V2Dt and the 
conventional WL result. 
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After restoring the scale dependence of the r's by R{t) 
we have, 



a{L) = (To 
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Evidently for R{Tie) ^ -R(t), cr{L) negatively diverges 
in d = 1,2, indicating that the back-scattering contribu- 
tion no longer has a perturbative effect on the a and such 
states are strongly localized. If one performs an expan- 
sion of the scaling function P{g) = ding / dlnR{Tie) near 
e = d-2 for e < 1, what follows is = d-2-l/q + . . . 
or the celebrated scaling theory of localization 21]. 

In general, r^e is a functional of R(t) and depends on 
the power spectrum of the environmental fluctuations. 
Experimentally, the result above can be detected in the 
limit of low temperatures T — ^ 0, given that r^e is ex- 
pected to vanish as a power law in T, or r^e (x T^. Con- 
sequently, in a Fermi-liquid a carefull determination of p 
is especially required in d = 2 where the obvious logarith- 
mic correction must be distinguished from the additional 
logarithmic correction in T generated by the exchange 
diagrams of the disordered-interacting problemQ- This 
can be accomplished by applying a small normal mag- 
netic field to suppress WL effects Q. Nevertheless, we 
have provided a self-contained description of the neces- 
sary mathematical formalism, indispensable for the rest 
of the paper. 

Interactions and Reaction- Diffusion- As we have 
shown ordinary diffusion with R{t) t^^"^ produces the 
WL corrections (|0.5p . We strongly emphasize that an al- 
ternative way of modeling a different groundstate can 
be accomplished by including a non-trivial ip{CD,r,t) 
in Eq. ()0.2p with possible non-linear and non-equilibrium 
terms. The generalized Fisher's equation which con- 
tains tp{CD,r,t) cx C_d(1 — Cd) is the simplest case of 
reaction-diffusion model, introducing local energy penal- 
ties or self interaction for Co coupled with competing 
growth terms, often useful in modeling problems related 
to quantitative biology and physical chemistry in addi- 
tion to non-equilibrium electronic systems[22]. In order 
to get an accurate sense of the scale dependent 5a and 
metal-insulator criticality, one needs to determine spe- 
cific forms of R{t) in arbitrary dimensions, which can be 
difhcult and is itself an active area of research. 

An approximate asymptotic solution of the Fisher's 
equation in arbitrary dimensions for a wide range of 
boundary conditions was performed by S. Puri, K.R. El- 
der and R.C. Desai23|, who determined the following 
expression for the spreading of the domain. 



interest, R{t) ^ t, and the 5a via (|0.4[) become. 
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for a fixed value of = 2, which is the velocity of 
the traveling wave solution of Cd- RF{d,t) takes on 
real values beyond the branch point of the squareroot at 
to = -{d/2)WL [-l/(2c;7r)], where Wl is the Lambert-W 
function. At very small values of t, R{t) is dominated by 
the logarithm, however at t » to, which is the limit of 
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- HR{t^)) d = 1 



— ^cx<( -2ft(i?(T^))-i d = 2 (0.8) 
-3ft2(i?(T^))-2 d^3 



which is clearly non-singular for large Tip and contains 
non-localized states in d = 2 for the weak disorder limit. 
This result is rather unexpected, given that we did not 
introduce any additional scaling couplings, rather a sim- 
ple unweighted quadratic term in ^/^(C_D,r, i) ^ {CdY is 
sufficient to circumvent localization. 

Interactions, Flory Scaling and Non Markovian Trans- 
port- Another unconventional groundstate having mem- 
ory dependent effects with generalized solutions of R{t) 
in arbitrary dimensions can be constructed by enforcing 
self avoidance onto the random scattering trajectories. 
Self avoidance will non-trivially introduce long range cor- 
relations into the impurity scattered trajectories of the 
electrons in contrast to directly considering correlations 
in the random potential itself. Although there might be 
a relationship between the two types of correlated ran- 
domness, here we are motivated to include repulsion be- 
tween the trajectories in order to induce the site or path 
of a scattered electron to have a dependence on its pre- 
vious scattering history [23|. Physically, a local repulsive 
term can be introduced into a random walk such that the 
walk never crosses itself and is free of intersecting loops 
as shown in FiglIJa). Such a model was pioneered by the 
Nobel prize winning chemist Paul Flory, who developed 
the famous scaling relation of a random walk with self 
avoidance [25|, 
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where the exponent associated with the spreading of the 
electron density will have an explicit dependence on the 
dimensions of the system. If we introduce Eq. (|0.9p into 
(lO^ we have 
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which is plotted in Figl2] and after restoring the scatter- 
ing length dependence for each value of d, our final result 
for the scale dependent a is given by 



ln(i?(l,v)/i?(l,r)) 



1 



cx <( 2(i?-2/3(2,T^) - i?-2/3(2,T)) d = 2 
f (i?r4/3(3,r^) - i?-4/3(3,r)) d = 3 

(0.11) 

showing exclusive localization only in d = 1, and the pos- 
sibility for a mobility edge in 2D now emerges near the 
band center. Although a 2DMIT is expected at some 
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FIG. 2. (Color Online) Conductivity corrections in a non- 
Markovian transport process arising from self-avoidant diflFu- 
sive scattering. Notice that the logarithmic behavior in d = 1 
is a clear indication of the lower critical dimension of delo- 
calization for which the perturbative /3 function in a 2 + e 
expansion fails to accurately capture the transition. 



critical value of the ratio of randoni potential width over 
the hopping strength, it is necessary for accurate numer- 
ical simulations using the finite size scaling method to 
confirm these aspects of the transition, and determine 
the critical exponents. 

Discussion and concluding remarks- For both cases 
studied here, we have observed a failure of localization in 
d — 2 simply by restricting the back-scattering by means 
of a modified diffusion process and determining its per- 
turbative effect on the transport. These two new types 
interactions in disordered electron systems may refiect an 
alternative ground state to the conventional Fermi liquid. 
First we have calculated the transport corrections in a 
Fisher's reaction- diffusion model, where the interactions 
are present as nonlinearities in the diffusion (or self in- 
teraction of the probability density) rather than directly 
taken as a density-density (or four fermion) interaction 
in the original Hamiltonian. 

Thus we have chosen as an alternative to a many-body 
picture, a one particle system with modified kinetics, 
which refiects a possibly different fixed point of the elec- 
tron liquid. Furthermore, in the Fisher's equation there 
is a local repulsion of the probability density in addition 
to non-equilibrium growth which may be an emergent 
manifestation of the strongly correlated electron liquid, 
as both theoretical and experimental evidence suggests 
that there may be density inhomogeneities, negative com- 
pressibility states 1^ 16 1 near the 2DMIT and the pre- 
dicted nematic phases that are intermediate between the 
Fermi-liquid and the Wigner crystal[26']. Hence, experi- 
mentally, by tuning the density of a 2DEG in a disordered 
SI-MOSFET heterostructure in the vicinity of « 10 
near the 2DMIT critical point we expect the negative 
compressibility conditions modeled by the phenomeno- 



logical considerations of a reaction-diffusion transport 
model to hold accurately. 

Regarding non-Markovian transport, recent cr(T) mea- 
surements in disorder tuned magnetic films have detected 
regimes with pure power-law like corrections to g(D(27|. 
Additionally ferromagnetic semiconducting systems are 
known to contain non-Markovian transport and spin re- 
laxation mechanisms for which anomalous temperature 
dependent corrections have been experimentally observed 
in d = 1,2(2^ [2^. Hence, materials such as ferromag- 
netic (Ga,Mn)As are likely candidates for observing the 
effects presented here. 

In conclusion, we have identified several unconven- 
tional modes of diffusive transport that can circumvent 
localization and it is shown here for the first time that 
non-equilibrium conditions can destroy the localization 
corrections in 2D. We attempted to develop a discus- 
sion of the single parameter scaling of f3{g) which should 
depend on the special symmetries of the Hamiltonian. 
Hence, as a compliment to our semiclassical scaling the- 
ory it is essential that several possible microscopic Hamil- 
tonians and stable fixed points that refiect the environ- 
mental conditions for these forms of restricted backscat- 
tering can be properly reconstructed and those relevant 
symmetries are identified in addition to careful numeri- 
cal simulations which are beyond the valid perturbative 
regime of our analysis. 
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